On the long range correlations of thermodynamic systems out of equilibrium 
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Experiments show that macroscopic systems in a stationary nonequilibrium state exhibit long 
range correlations of the local thermodynamic variables. In previous papers we proposed a Hamilton- 
Jacobi equation for the nonequilibrium free energy as a basic principle of nonequilibrium thermo- 
dynamics. We show here how an equation for the two point correlations can be derived from the 
Hamilton- Jacobi equation for arbitrary transport coefficients for dynamics with both external fields 
and boundary reservoirs. In contrast with fluctuating hydrodynamics, this approach can be used 
to derive equations for correlations of any order. Generically, the solutions of the equation for the 
correlation functions are non-trivial and show that long range correlations are indeed a common 
feature of nonequilibrium systems. Finally, we establish a criterion to determine whether the lo- 
cal thermodynamic variables are positively or negatively correlated in terms of properties of the 
transport coefficients. 
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The basic idea in the construction of nonequilibrium 
thermodynamics is that of local equilibrium. This means 
that on the macroscopic scale it is possible to define lo- 
cal thermodynamic variables like density, temperature, 
or chemical potentials, which vary smoothly. Microscop- 
ically this implies that the system reaches local equilib- 
rium in a time which is short compared to the times typ- 
ical of macroscopic evolution as described, for example, 
by hydrodynamic equations. Even if the local behavior 
of a nonequilibrium system is well described by a Gibbs 
state, if we probe the system at distances comparable 
to its macroscopic size, relevant differences with respect 
to equilibrium do appear. It is indeed an established 
fact that macroscopic systems in a stationary nonequi- 
librium state generically exhibit long range correlations 
of the local thermodynamic variables. For a review of 
experimental results see The theories of correlations 
in thermodynamic systems out of equilibrium that can 
be found in the literature are mainly based on the so 
called fluctuating hydrodynamics 0, S 0| , that is on a 
linearization around a steady state of the hydrodynamic 
equations perturbed by random currents. For diluted 
systems this problem has also been approached within 
kinetic theory Q. 

It has been shown for several stochastic dynamics that 
the nonequilibrium free energy, which is the generat- 
ing functional of the correlation functions, satisfies a 
Hamilton- Jacobi equation 0, Q. The basic feature of 
this equation is that it does not depend on the details of 
the microscopic dynamics but only on the macroscopic 
transport coefficients, i.e. the diffusion coefficient and 
the mobility. It can thus be considered as a general prin- 
ciple of nonequilibrium thermodynamics. In this paper 



we show how the Hamilton-Jacobi equation can be used 
to study correlation functions. In contrast with fluctuat- 
ing hydrodynamics, this approach can be used to derive 
equations for correlations of any order. For simplicity, 
we consider diffusive systems in arbitrary dimension de- 
scribed by a single thermodynamical variable with both 
external fields and boundary reservoirs. We show that 
the two point function satisfies a linear partial differen- 
tial equation which depends, as in fluctuating hydrody- 
namic, on the transport coefficients evaluated in the sta- 
tionary state. These coefficients can be computed from 
the microscopic dynamics or can be obtained from exper- 
iments. Generically, the solutions of the equation for the 
correlation functions are non-trivial and show that long 
range correlations are a common feature of nonequilib- 
rium systems. We remark that the long range correla- 
tions arise from nonlinear effects in the hydrodynamic 
equation. One of the most interesting consequences of 
the equation for the correlations is that it is possible to 
establish whether the fluctuations are positively or neg- 
atively correlated in terms of structural properties of the 
transport coefficients. 

Systems for which the above scheme can be applied are 
characterized by a separation of scales both in space and 
time. We analyze these systems in the limit in which this 
separation becomes sharp. In the following we let e be 
the ratio between microscopic scale, say of order of the 
intermolccular distance, and the linear size of the system 
Lo- Let us consider a macroscopic system in a domain A 
in space dimension d described in terms of a single ther- 
modynamic observable /?, that we think of as the density 
of particles. We denote by p(x) the local density at the 
macroscopic point x, obtained as the density of parti- 
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cles in a volume of side length £, with eLq <C t <C Lo- 
We then introduce the corresponding ensemble P e as the 
distribution of such local density. This distribution is 
inherited from the microscopic ensemble. 

According to the Boltzmann-Einstein theory of equi- 
librium thermodynamic fluctuations Q the relationship 
between the free energy F and the probability distribu- 
tion P e is the following 



P e {p{x),x^k)ote- F[p)/£d 



(1) 



where we absorbed the temperature T, fixed in the one 
phase region, in the definition of F . In this equilibrium 
setting the free energy F is a local functional, i.e. it can 
be written as 



Hp) 



dxf(p(x)) 



(2) 



where f(p) is the free energy per unit volume. In par- 
ticular, there are no correlations, i.e fluctuations of p(x) 
and p(y), for x ^ y, are independent. 

For systems out of equilibrium, e.g. when they are in 
contact with reservoirs and/or under the action of an 
external field E, we let P £ be the corresponding station- 
ary ensemble. For such systems we define the associated 
noncquilibrium thermodynamic functional F(p) by (fT]). 
That this is the appropriate definition is supported by 
the following facts discussed in 0, S|- First, the sta- 
tionary state corresponds to the absolute minimum of 
F, thus generalizing the maximum entropy principle of 
equilibrium statistical mechanics. Moreover, the deriva- 
tive of F with respect to p provides the correct definition 
of thermodynamic force, that is the force responsible for 
dissipation, as shown in Q. 

In order to present the Hamilton- Jacobi equation for 
the free energy, we introduce the dynamical behavior of 
the system. The macroscopic evolution of the density is 
described by a (in general nonlinear) diffusion equation 
with a transport term corresponding to the external field 
E, namely 

d t p(t,x) = di -D ij (p)d j p(t,x) - X ij{p)Ej{x) (3) 

where we use Einstein's convention of summing over re- 
peated indices and di stands for d/dxi, i = 1, • • • , d. In 
([3]) D is the diffusion coefficient and \ is the mobility. 
This equation has to be supplemented by the appropri- 
ate boundary conditions. For noncquilibrium systems in 
contact with particle reservoirs, this amounts to fix the 
value of p at the boundary of A. We denote by p(x) 
the stationary density profile, i.e. the unique station- 
ary solution of ([3]). This is the density profile in which 
the noncquilibrium free energy F attains its minimum. 
Equation ([3]) identifies the current J of particles flowing 
though the system in the stationary profile as 



Ji(x) 



- Dijipjdjpix) - Xij(p)Ej(x) 



(4) 



The equation for p then reads diJi(x) — 0, i.e. J is di- 
ver genceless. 

As shown in 0, @] , the nonequilibrium free energy F is 
the maximal solution of the Hamilton-Jacobi equation 



SF 1 
Tp' 2 



Xij{p)d- 
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(p)d j p + Xij (p)E j )=0 



(5) 

where (/,<?) = J A dx f(x)g(x) and F satisfies the bound- 
ary condition SF/6p\ g ^ = 0. The derivation of this 
equation is based on a dynamical argument. One first 
generalizes the Boltzmann-Einstein formula (TTJ) to space- 
time trajectories obtaining the asymptotic probability of 
deviation from solutions of the hydrodynamic equation 
The functional F can be identified as follows: if at 
t = — oo the system is in the stationary profile p then 
the probability of observing the profile p at t = is pro- 
portional to exp{— e~ d F(p)}. The free energy F solves 
a variational problem whose associated Hamilton-Jacobi 
equation is ((5|). The validity of ([5|) has been established 
for stochastic lattice gases for which the local equilibrium 
can be proven. 

In the case of equilibrium systems, in which E = and 
p\dA is constant, the profile p is constant; it is then simple 
to check that the solution of (0 has the form |J2J) where 
the free energy per unit volume f(p) satisfies the Ein- 
stein relation Dij(p) = Xij(p)f"(p)t normalized so that 
f'(p) = f(p) = 0- For nonequilibrium systems the solu- 
tions of ([5]) cannot, in general, be obtained in a closed 
form. For special choices of the transport coefficients, 
which correspond to well studied stochastic dynamics, 
the solution of ([5]) can be expressed in terms a non linear 
boundary value problem, see [1, [l^, 11, 12 1. Contrary 



to equilibrium, a general feature of the noncquilibrium 
free energy F is that it is not a local functional. This im- 
plies that macroscopic density correlations do gcncrically 
appear. 

The density-density correlation is defined as 



C(x,y) 



(p(x)~p(x))(p(y)-p(y))dP £ (6) 



Note that the scaling e~ d corresponds to the Gaussian 
fluctuations. This function is related to the non equilib- 
rium free energy F by 



C-^x^y) 



S 2 F(p) 



5p{x)5p(y) 



p=p 



(7) 



so that F(p) = (1/2) ((p - p), C-\p - p)> + o((p ~ p) 2 ). 
We next introduce the function B(x,y) by 

C(x,y) = C eq {p{x))6(x-y) + B(x,y) , x,y e A (8) 

where, as / is the equilibrium free energy per unit vol- 
ume, C eq {p) = /"(p) _1 gives the (local) equilibrium vari- 
ance. Notice that, while D and x arc d x d matrices, C eq , 
C and B are multiples of the identity. Since the density 



3 



fluctuations at the boundary are determined only by the 
reservoirs, we have that B(x,y) = when either x or y 
is at the boundary of A. 

In order to write an equation for B, introduce the el- 
liptic differential operator L as 

L=\ Dij (p(x)) 6^ + X 'ij (p(x))Ej (x) d t (9) 



We now consider some onc-dimensional systems and 
choose A = (0, 1). We let p(0) = p < p\ — p(l) be the 
boundary conditions imposed by the reservoirs. The one- 
dimensional exclusion process has < p < 1, D(p) = 1, 
x(p) = p{l — p). For E = equation (flU)) implies that 
B is proportional to the Green function of the Dirichlct 
Laplacian in (0, 1), namely 



and£ = L©L, i.e. C satisfies £<p(x)ip(y) = i/>(y)L<p(x)+ B(x, y) = -(pi - p ) 2 x(l -y) , < x < y < 1 (13) 
(p(x)Lip(y)- Let finally L* be the adjoint of L with Dirich- 
let boundary conditions on A, i.e. 



^VO) = -didj [Dij (p(x))<p(x)] 

-d^ip^mixMx)] 

and = © . As shown below, from the Hamilton- 
Jacobi equation (JSJ) it follows that B satisfies 



tfB(x,y) = -h(x)5(x-y) 



(10) 



where 



h(x) = -di[y!y{p(x))Dj£{p(x)) J k (x)] (11) 



The choice of the Dirichlct boundary condition for L cor- 
responds to the vanishing of B at the boundary. 

We draw some important consequences of equation 
(flO]) . Since D is positive definite, C is an elliptic op- 
erator in A x A with Dirichlct boundary condition. Let 
G(x,y; ,x',y') be its Green function, i.e. the solution of 
— CG(x,y;x' ,y') = S(x — x')S(y — y'). Then the solution 
of HDD is 



B(x,y) = / dz G(z, z; x,y) h(z) 



(12) 



Since G > 0, wc conclude that if h > 0, respectively 
h < 0, fluctuations of the density arc positively corre- 
lated, i.e. B > 0, respectively negatively correlated, i.e. 
B < 0. In absence of external field, i.e. for E = 0, wc 
have, from ^ and (flT|) . that h(x) = ^didjXij(p(x)). If 
we further assume that x is a multiple of the identity, 
we get that B > 0, respectively B < 0, if \{p(x)) is a 
subharmonic, respectively superharmonic, function of x, 
namely Ax(p(x)) is positive, respectively negative. 

For equilibrium systems the current J in the station- 
ary profile vanishes, hence h = so that C = C eq . This 
conclusion is also true for systems with periodic bound- 
ary conditions and constant external field. In such case 
p and J are constant so that h vanishes. More gener- 
ally, in (l3l ] it is shown that, in the periodic case, the 
nonequilibrium free energy F is the same as the equilib- 
rium one. For isotropic systems D and x are multiples 
of the identity, in such a case h vanishes if and only if 
&i [x' {p{ x ))D~ 1 {p{x))\ Ji{x) = for any x in A. For in- 
stance, in the zero range process x is an increasing func- 
tion and D = x', so that h vanishes and the model does 
not exhibit long range correlations. 



This result was first derived in Q. Notice that, in agree- 
ment with the above discussion, since x is concave and 
p' is constant, B < 0. For constant external field E the 
solution of (flQ|) is given by 



B{x,y) = 2Jp'{x)p\y) 



J*dup'(u) 1 Jydup'(u) 1 
fodup'iu)- 1 



(14) 



for < x < y < 1. This formula for the correlation 
function has been derived in ljj. It can be shown that 
p is increasing. Therefore the correlations B(x,y) have 
the same sign as J. We finally discuss the correlations 
in the Kipnis-Marchioro-Presutti model [HI, [HI]. This 
model describes a chain of one-dimensional harmonic os- 
cillators which are mechanically uncoupled and interact 
by exchanging stochastically the energy with the neigh- 
boring sites. Accordingly, the thermodynamic variable 
p is the energy density. For this model D{p) = 1 and 
x(p) = p 2 ■ The solution of (flQ|) is then the same as 
the one for the exclusion model, i.e. (fl"3"|) . but with the 
opposite sign. For this model the energy fluctuations 
are therefore positively correlated, in agreement with the 
above discussion and the fact that x is a convex function. 
Since B is proportional to Green function of the Dirichlct 
Laplacian, we also have that it is a positive definite op- 
erator. This means that, in the quadratic approximation 
near p, the nonequilibrium free energy is smaller than the 
local equilibrium one. The opposite behavior takes place 
in the exclusion process. For these models this result has 
been proven also for large fluctuations, that is beyond 
the quadratic approximation 0, 1^, 11|. 

We here show that the above results can be seen as 
special cases of a class of one-dimensional systems for 
which (|10p can be solved explicitly. Assume that 



2x'(p{x))E(x)h{x) = d x [D(p(x))h(x)\ (15) 

holds for any x in (0, 1). Then it is simple to check that 
the operator L = —2h(x)~ 1 is self-adjoint and B is its 
Green function. By standard Sturm Liouville theory we 
then get 



B{x,y) = ^a 1 (x)a 2 {y), 



0<x<y<l (16) 



where a^, i = 1,2, solves the Cauchy problem Lai = 0, 
with ai(0) = 0, ai(0) = 1 and a 2 (l) = 0, a' 2 (l) = 1. 
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Finally K = -[D(p(x))/h(x)]W(ai,a 2 ), where W is the 
Wronskian, is constant in x. Both for the exclusion and 
the KMP processes (fT5)) holds for any constant external 
field E. Simple computations then yield (|13[) and (fT4"|) 
as special cases of ([TJJ]). Moreover, from (flTJl) we also get 
the correlation function B(x,y) for the KMP model with 
constant external field, which is the same as ([T4]) with 
opposite sign. 

To conclude we show how equation (jlOj) can be derived 
from the Hamilton-Jacobi equation for the nonequilib- 
rium free energy. To derive the equation for the correla- 
tions, it is convenient to introduce G(A) as the Legcndre 
transform of F(p), i.e. G(A) = sup p { (A, p) — F(p)}. The 
function X(x), x in A, can be interpreted as the variation 
of the chemical potential, in particular A = corresponds 
to the stationary profile p. By Legendre duality, equation 
([5]) is then equivalent to 



that is (TlPl). 



d t X,l D J^)d. 



5G 



(8G 







for any A which vanishes at the boundary of A. In 
the quadratic approximation we have G(A) = (A, p) + 
(1/2) (A, CX) + o(A 2 ). Since p is the stationary solution 
of ([3]), the first order term in A vanishes in the above 
equation, so that G solves 

a, a, \xiM^ - \d 3 (D t3 {p)cx) +x! ij {p)E j cx) = o 

By plugging (J8]) in the above equation, recalling ([4]) and 

D Ik(P)Xk 3 {p) = C eq (p)Sij, we get 



9iX , --d J {D lJ (p)BX)+ X ' lJ (p)E J BX 



diX^iMDT^p)!^ 



Integrating by parts we then get 



A , \didi (Dij (p)BX) - ft (x'a (P)E 3 BX) 



diix'ijiP^IkJk) 



Recalling ([9]) and (JTTJ , the above equation can be written 
in an operator notation as 

L^B + BL = -h 

Here h denotes a multiplication operator and we used 
that B{x,y) = B(y,x). Denote by L\, respectively Lj, 
the operator L' acting on functions of x, respectively of y. 
From the previous equation we get that for any function 
ip vanishing on 9A we have 

/ dy[LlB{x,y)tp{y) + LlB{x,y)tp{y)] = -h(x)cp(x) 



We conclude with some comments on the applicabil- 
ity of equation (10). The Hamilton-Jacobi equation from 
which it is derived was originally introduced in the study 
of microscopic models whose macroscopic behavior is of 
diffusive type and there are no normal modes like sound 
waves. This case covers for example nonequilibrium so- 
lutions of different chemical species but no chemical re- 
actions. The Hamilton-Jacobi equation is strictly con- 
nected with the dynamic generalization of the Einstein 
formula for static fluctuations obtained in 0, Q and 
it has been extended to systems without conservation 
laws [TH, [3] where the hydrodynamic equations are of 
reaction-diffusion type, as it happens in presence of chem- 
ical reactions. It can be generalized to situations with 
several conservation laws where one of the hydrodynamic 
equations is Navier-Stokes. We finally emphasize that the 
approach here presented, based on an exact equation for 
the nonequilibrium free energy which takes into account 
the nonlinear effects in the systems, can be applied to 
derive equations for higher order correlations. 
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